The Exponential Matrix 



The work in the preceding note with fundamental matrices was valid 
for any linear homogeneous square system of ODE's, 

x' = A{t)x. 

However, if the system has constant coefficients, i.e., the matrix A is a con- 
stant matrix, the results are usually expressed by using the exponential ma- 
trix, which we now define. 

Recall that if x is any real number, then 

= 1 + ^ + 1^+... + ^ + ... . (1) 
2! n\ 

Definition 3 Given an n x n constant matrix A, the exponential matrix e"^ 
is the n X n matrix defined by 

= j + A + ^ + + ^ + . (2) 

Each term on the right side of (2) is an n x n matrix adding 
up the zy'-th entry of each of these matrices gives you an in- 
finite series whose sum is the ij-th entry of e"^ . (The series 
always converges.) 

In the applications, an independent variable t is usually included: 

e^^ = J + A f + ^ + . . . + A" + . . . . (3) 
2! n! 

This is not a new definition, it's just (2) above applied to the matrix A Mn 
which every element of A has been multiplied by t, since for example 

(Atf = At -At = A-A-t^ = A^t^. 



Try out (2) and (3) on these two examples (the second is very easy, since 
it is not an infinite series). 

Example 3A. Let A = ^ ^ ^ ^ . Show: e"^ = ^ ' 
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Example 3B. Let A = ( ^ Q ) ' show: = f J j ) and 
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What's the point of the exponential matrix? The answer is given by the 
theorem below, which says that the exponential matrix provides a royal 
road to the solution of a square system with constant coefficients: no eigen- 
vectors, no eigenvalues, you just write down the answer! 

Theorems Let A be a square constant matrix. Then 

(1) (a) e^^ = <i>o{t), the normalized fundamental matrix at 0; 

(2) (b) the unique solution to the IVP V = Ax, x(0) = xq is x = 

e^^XQ. 

Proof. Recall that in the previous note we saw that if Oq (f) is the normal- 
ized fundamental matrix then 

The solution to the IVP : x' = A{t) x, x(0) = xq is x{t) = ^o{t)xo. 

(4) 

Statement (2) follows immediately from (1), in view of (4). 

We prove (1) is true by using the fact that if = then the normalized 
fundamental matrix has <t>(0) = I. Letting <t> = e"^*, we must show = 
A<I> and<I>(0) = 1. 

The second of these follows from substituting t = into the infinite 
series definition (3) for e'^^ 

To show O' = AO, we assume that we can differentiate the series (3) 
term-by-term; then we have for the individual terms 

A« ^ = A" ■ 



since A" is a constant matrix. Differentiating (3) term-by-term then gives 

= j-^e^' = A + A^t + ... + A" ^ 
Ae^* = AO 



dt - dt'- - ^-L -r I -r . . . -r („_i^, -r . . . 

— A a^t — ' ^ 



Calculation of e'^K 

The main use of the exponential matrix is in Theorem 3 — writing down 
explicitly the solution to an IVP. If e^'^ has to be calculated for a specific 
system, several techniques are available. 
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a) In simple cases, it can be calculated directly as an infinite series of 
matrices. 

b) It can always be calculated, according to Theorem 3, as the normal- 
ized fundamental matrix $o(0' using (11): Oo(i) = 

c) A third technique uses the exponential law 



(6) 



valid if BC = CB. 
To use it, one looks for constant matrices B and C such that 

A = B + C, BC = CB, e^* and e^^ are computable; (7) 



then 

Example 3C. Let A = 

using e^K 

Solution. We set B = 

satisfied, and 



e^' = 

2 1 
2 

2 
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(8) 



. Solve x' = Ax, x(0) = ^ ^ 

then (7) is 



and C = 
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by (8) and Examples 3A and 3B. Therefore, by Theorem 3 (2), we get 



X = e^^xo 
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